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I.  IHTRODUCTION 

General 

Deposition  of  aerosol  particles  on  various  obstacles  is  of  extreme 
importance  in  many  natural,  technological,  and  scientific  systems.  For 
example,  performance  of  industrial  installations  for  abatement  of  parti- 
culate pollutants  (spray  scrubbers,  cyclones,  electrical  precipitators, 
and  fibrous  filters),  air  pollution  field  samplers  (impactors,  small-pore 
filters),  and  the  human  lung, depends  on  the  deposition  characteristics  of 
the  particles. 

Analyzing  the  behaviour  of  aerosols  in  the  collecting  systems,  one 
is  confronted  with  the  significant  difference  between  spherical  and  non- 
spherical  particles.  The  difference  is  mainly  due  to  the  poor,  or  even 
lack  of  symmetry  of  the  latter  (Gallily,  1974),  which  makes  the  under- 
standing of  their  aerodynamic  behaviour  very  difficult. 

To  achieve  a workable  knowledge  of  the  aerodynamic  behaviour  and 
the  deposition  rate  of  the  particles  on  obstacles,  one  usually  follows  two 
lines: 

a.  Theoretical  solution  of  the  equations  of  motion  for  typical, 
illuminative  cases; 

b.  Experimental  studies  of  these  cases. 

Obviously,  the  lines  have  to  complement  each  other,  thus  establishing 
confidence  in  the  method  of  problem  treatment. 


Equations  of  motion 


Depending  on  the  size  of  the  particles  and  the  nature  of  the 
flow  which  carries  them,  one  can  deal  with  the  motion  of  the  foruer  from 
a deterministic  or  a stochastic  point  of  view. 

In  the  first,  one  has  to  solve  the  equations  of  translation  and 

rotation,* 

™ = I rt(i  + Fr  Hi] 

i 

and 


Zl 


n,t,  l 


+ L 


cn 


whereas  in  the  second,  one  has  to  deal  with  the  equation  of  convective 
diffusion 


dcx/ bt  + C V + B-Fe  ) • VC,  r cl'iv  ( \*T  B ‘ VcJ  C 

which  are  subjected  to  initial  and  boundary  conditions  typical  of  the 
physical  situation. 

The  main  difficulty  in  the  solution  of  equations  [ll-C?]  lies  with  the 
evaluation  of  the  fluiddynamic  force,  torque, and  mobility  tensor,  formally 
defined  for  a continuous  medium  by 

f'»Ks-2r  ’ m 

g 

= J >XC«»S-TT)  , C5J 

S 


For  symbols  see  GLOSSARY 


™ 


and  (Happel  and  Brenner,  1965a) 

B = C-S  ciS-P]  *■  U3 

/y  r w 

where  TT  and  P are  respectively  the  stress  diadic  and  triadic  related  to  the 
field  of  flow.  Thus,  the  knowledge  of  this  field  is  a prerequisite  to  the 
solution  of  the  equations  of  motion. 


Scientific  background  and  previous  studies 


A.  Mobility.  Forces,  and  Torques 

The  aerodynamic  force  and  torque  acting  on  particles  embedded  within 
an  arbitrary  field  of  flow  near  a collecting  wall  has  been  theoretically 
calculated  and  experimentally  determined  for  a few  types  of  cases. 

There  is  a limited  number  of  studies  concerning  a motion  of  a spherical 
solid  particle  (Brenner,  1965;  Brenner,  Goldman,  and  Cox,  1967;  Gallily  and 
Mahrer,  1973),  and  even  less  concerning  a motion  of  spheroids  (Wakiya, 

1959)  or  cylinders  (Be  Mestre  and  Russel,  1975),  which  take  care  of  the 
effect  of  the  wall.  The  overwhelming  majority  of  investigations  were 
devoted  to  the  study  of  the  fluiddynamic  force  and  torque  which  operate 
on  isolated  particles,  assuming  that  the  interaction  between  the  latter 
and  the  wall  can  be  neglected.  In  these  studies,  the  field  of  flow  was 
characterized  by  the  creeping  motion  conditions 

«•«  = LW|  /t  « 1 CV 


= Le  I “I  /s>  «i 
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where  f is  a typical  length  of  the  particle,  u)  is  the  rotational 
velocity,  and  is  the  kinematic  viscosity. 


1.  Motion  of  isolated  particles  in  a continuous-regime 
a.  Ellipsoids 

Considering  the  case  of  relatively  large  particles  for  which  Knudsen 
number  Kn  is  Kn  (=  A / £ ) <£<  i , one  has  for  a steady  motion  of  an 
isolated  ellipsoid  within  an  arbitrary  flow  v (Brenner,  1964) 

= r **  t + • (j^v)  - •••> 

<V  ^ 3 , D 1 ' t 

~ 1 ral 


l„,p  = f*  | ^ (Daa  x v)„  + ~j~  (c"dxv), 

+ -IvA.-  (x>fc  a xv)M  + •••>-  (a*.,  + a^2)  cO]  [103 

where  K , the  translation  tensor,  is 
* 

K = l<B  TV  ft  C.  — ii + ~ ~ \ rill 

* ,1'  **-«?*•  *.+&.*  * 


Q.  is 


a = -!4e«A s( 4-A  , . i i.  , it 

* 3 1 J c ]’ 

x*" 

O0 

r-  - «&./  ■*. 

"D  “ ^ / £>*'*  ■+•  6 , dZ/  d^,A  + cA  £>*/  ; 


are  the  ellipsoid's  semi- 


i afS/ix'-j  «?>/»»'.  c,,  6-,,  c, 

axes;  and  i ■ 1,2,3  stands  respectively  for  a,b,c,  in  a cyclic  order. 

For  a non-steady  motion  of  an  isolated  particle  in  an  arbitrary  non- 
steady field,  one  has  to  take  into  account  an  "added  mass"  and  "added 
moment  of  inertia"  (Lamb,  1932),  an  "added  buoyancy",  and  relaxation  effects 
resulting  from  the  changing  relative  velocity  between  the  particle  and  the 
fluid  (Morrison,  1974),  The  "added  mass"  and  "added  moment  of  inertia"  which 
have  to  be  combined  with  the  mass  and  moment  of  inertia  of  the  particle, 
respectively,  are  proportional  to  the  density  of  the  displaced  fluid,  the 
proportionality  factor  being  of  the  order  of  unity.  So,  for  aerosol  systems 
they  can  be  neglected.  On  the  other  hand,  the  so  called  Basset-term, 
related  to  the  fluid-relaxation  effects  has  yet  to  be  checked  for  non- 
spherical  particles. 

In  the  case  of  an  isolated  solid  sphere,  for  example,  the  most  general 
expression  of  the  fluiddynamic  resistance  F>.  is  (Morrison,  1974e) 


_ Gv 


ra(u:-v;)-  ^ c 

— w (t  — t^) 


where  aois  the  radius,  m'  is  the  mass  of  the  displaced  fluid,  and  J>/srt  is 
the  substantial  derivative  operator.  Thus,  the  introduction  of  [l2]  into 
the  equation  for  the  creeping  translational  motion  of  even  a solid  sphere, 


would  turn  the  last  equation  into  an  integro-differential  one.  This 


introduction  has  been  never  done  due  to  the  enormous  complications  it 
brings  to  the  equations  of  motion  [l 2-12]  for  even  the  simple  spherical 
particle* 

Here  it  should  be  remarked  about  the  suitability  of  an  ellipsoid  to  simulate 
the  shape  of  real  particles.  Approximating  the  shape  of  the  latter  by  that 
of  regular  bodies,  one  can  distinguish  among  categories,  as  given  in  Table  1 


Shape 


Spherical 

Irregular 

Cubical 

Flakes 

Fibrous 

Condensation 


Percent  by  Weight 
in  sample 


Range 


Average 


0-20 

10 

10-90 

40 

. 

0-10 

5 

5-35 

1 

10 

0-40 

15 

Table  Is  Airborne  Dust  Particle  Shapes  (stern,  1968) 
(*  Also,  asbestos  fibers  and  cloud  ice  crystals). 


Particles 
(larger  than 
0.1  m) 


Smoke,  pollen 
fly-  ash 

Mineral 

Cinder 

Mineral 

epidermis 

Lint  plant 
fiber* 

Carbon,  smoke 


* 
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In  this  category  system,  a cylindrically  shaped  particle  (fiber,  needle) 
may  be  looked  upon  as  a prolate  rotational  ellipsoid  of  a very  large 
polar  to  equatorial  axis  ratio  or  a finite  slender  body.  A discoid-type 
particle  (flake)  may  be  envisaged  as  an  oblate  rotational  ellipsoid  of  a 
very  small  axial  ratio,  and  a straight-chain  condensation  floe  of  spheri- 
cal primary  particles  may  be  roughly  considered  as  a rod. 

For  the  asymptotic  forms  of  the  cylinder  and  the  disk,  one  has  instead 
of  Eqs.  [9]  and  Lioj  somewhat  simpler  equations  relating  k and  Q.  to  the 
axial  ratio  of  the  body  (Brenner,  1 965a;Gallily,  1975). 

b.  Regular  bodies  with  planar  surface  sections 

There  is  no  theoretical  solution  for  the  fluiddynamic  resistance  force 
and  torque  acting  on  particles  having  planar  surfac-  sections  such  as  cubes, 
etc.  (Brenner,  1965h;  all  studies  of  these  particles  have  been  empirical 
(Heiss  and  Coull,  1952;  Pettyjohn  and  Christianson,  1946  ; Jayaweera,  1969; 
KajKawa,  1971;  Horwath,  1974;  Walkenhurst,  1976).  Heiss  and  Coull  ( 1 952 ) , 
for  example,  defined  an  empirical  factor  G relating  the  rate  of  fall  a non- 
spherical  particle  to  that  of  a volumetrically  equivalent  sphere  of  diameter 
, by 

Uoo  = (2  § 9 cC/lS  vO  G"  C 1 33 

for  a flatwise  fall  and  constant  circularity  0(4/0/ , where  oln  is  the  diameter 
of  a circle  equal  to  particle's  cross  section  perpendicular  to  its  motion, 
they  found 


= ^,0  CdS^V0<-)-^2S-y'/2c/a/o/'r,,  pq] 


and  for  an  edgewise  fall  and  a variable  circularity, 


in  which  S'  , the  sphericity  factor,  is  S'  = 77  /s  and  & is  the  surface  of 
the  particle. 

The  relationship  between  G and  Kj;  is  given  by 


X il  = (3  t r/a^,  s')  cU 


C 163 


It  should  be  remarked  that  no  empirical  formula  is  available  for  the 
rotational  mobilities  of  the  discussed  bodies. 


2.  Motion  of  isolated  particles  in  the  molecular  regime 
The  mobility  of  aerosol  particles  in  creeping  motion  and  the  molecular 
regime  where  Kri » 1 has  been  calculated  for  a rigid  sphere  (Epstein,  1924), 
and  later  for  nonspherical  particles  of  a cylindrical,  discoidal,  spheroidal 
and  cubical  shape  (Dahneke , 1973a;  Schwartz,  1976).  Dahneke1 s general  equation 
for  the  resistance  force  in  a direction  characterized  by  the  cosines  o t1 
with  respect  to  the  velocity  uMis 

r*  = (?  %*/ Lc^)  ■*W4(‘rJ-r)'/2oiol' 

+ 0/2)1  PP'+O/O-f JT1  J sfS  tnj 


* 


A general  formula  only. 
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where  Lc  is  a characteristic  length,  = A/Lt  t an<*  f is  the  fraction 
of  gas  molecules  reflected  diffusely  with  a temperature  accommodation 
coefficient  of  1. 

3.  Motion  of  isolated  particles  in  the  continuum-slip  regime 
Aerosol  particles  in  many  natural  and  man-made  systems  have  sizes 
characterized  by  0 C Kn  <.  0.2;  however,  there  are  but  a few  studies  of 
mobility  of  nonspherical  particles  in  this  Knuds en-number  range. 

In  analogy  with  the  case  of  two  rigid  spheres  (Morrison,  1974b)  in  which 
the  slip  boundary  condition  on  their  surface  is  taken  to  be 

v*  “ ( */ r)  ^ His] 

where  C^sfe-'Q/f  and  <Jt  is  the  tangential  stress  on  the  surface,  one 
can  assume  that  a resistance  tensor  K calculated  for  a nonspherical  particle 

/v 

on  a continuous  fluid  assumption  can  be  adjusted  for  the  continuum-slip 
regime  by  multiplying  it  by  a factor  a which  depends  on  Kn  and  C^,  alone, 

• <*  CKw,  . Cl''] 

/V 

The  mobility  of  non  -spherical  particles  in  the  continuum-slip  regime  was 
calculated  oy  a semi-empirical  technique  (Dahneke,  1973).  In  this  technique 
one  finds  first  a sphere  of  radius  R having  the  same  ratio  of  resistances 
in  the  continuous  and  molecular  regime  as  that  of  the  particle  in  question; 
then,  assuming  that  sphere  to  have  also  the  same  ratio  of  the  continuous- 


regime  resistance  to  the  resistance  at  any  Knuasen  number  as  that  of  the 
nonspherical  particle,  one  deduces  from  the  equation 


- 10  - 


K / = Lc  / R 

and  the  typical  length  Lc  the  appropriate  ^^(sA/r)  and  calculates 
the  mobility  of  the  nonspherical  particle  by 

£>  — £>  co-rt.  C 1 + +■  £*"  ) , C20j 

In  the  experimental  side,  the  translational  mobility  in  the  continuum- 
slip  regime  was  studied  at  our  laboratories  (Gallily,  1975)  for  cylinders 
sedimenting  at  an  atmospheric  pressure. 


B. 


A directed  study  of  deposition  of  nonspherical  particles  on  various 
obstacles  has  not  been  performed  until  recently  when  Gallily  and  'coworkers 
(Gallily,  1975),  and  Gentry  and  Spurny  (1976)  have  started  investigations. 

The  rate  of  sedimentation,  however,  has  been  used  for  quite  a time  as 
a means  for  size  characterization  through  the  concept  of  the  Stokes  (Rs) 
or  the  Aerodynamic  ( Radius  (Fuchs,  1964;  St8ber,  1972;  and  others) 
respectively  defined  by 

r,  = (’M  £-.«!/ 2 1:3.1  ffyh 

or  cm 

R»  = ( 9 p IS”.«I  / 2|3  lsr)  ^ 


1 . Aerodynamic  radius  studies 

Studies  on  the  aerodynamic  radius  of  nonspherical  bodies  were  carried 
out  for  aggregates  of  spherical  latex  aerosol  particles  of  various  confi- 
gurations and  for  asbestos  fibers,  mainly  by  the  aid  of  the  Stttber  spectro- 
meter (Sttfber,  1972),  and  for  variously  shaped  irregular  aerosol  particles 
(Kotrappa,  1972),  all  in  an  atmospheric  pressure  and  Knudaen  numbers  of 
Kn  < 0.2,  approximately.  Also,  aerodynamic  radii  of  aerosol  latex  aggregates 
have  been  determined  at  reduced  pressures,  i.e.  higher  Knudsen  numbers 
(Hochrainer  and  HMnel,  1975). 

Experiments  on  the  rate  of  sedimentation  of  (model)  cylinders,  ellipsoids, 
cubes  etc.  in  liquids  at  creeping  flow  conditions,  have  been  carried  out 
lately  too  (Horvath,  1974;  Walkenhdrst,  1976), 

These  studies  were  concerned  with  isolated  particles  falling  far  removed  from 
walls;  no  study  was  dedicated  to  investigate  the  aerodynamic  effect  of  the 
wall.  The  result  for  aerosol  particle  systems,  where  obviously  ho  ideally 
regular  bodies  could  be  produced,  were  expressed  in  the  form  of  empirical 
equations  related,  in  the  case  of  linear  aggregates  and  fibers,  to  the 
theoretical  continuous-regime  equation  and  incorporating  a Cunningham- like 
correction  factor. 

2.  Diffusions 1 transport 

The  diffuaional  transport  of  nonspherical  particles  was  theoretically 
studied  by  Brenner  (1972)  by  treating  an  ensemble  rather  than  individual 


bodies, 
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3.  Stochastic  motion  of  single  particles 

The  sedimentation  in  still  air  of  nonspherical  particles  in  the 
micron-size  range  was  studied  theoretically  by  a Monte  Carlo  technique 
(Gallily,  1975;  Gallily  and  Cohen,  1976)  in  which  only  regular  shapes  such 
as  cylinders,  disks,  and  rotational  ellipsoids,  were  treated.  In  this  study 
the  probability  density  function  for  an  orientation  change  of  a single  sedimenting 
particle  was  taken  to  be  (Gans,  1928) 

P(X,+)  = f e. P„(X)  t25] 

n=o 

where  /=  cos  9 and  (X)  is  Legendre's  polynomial  of  degree  n. 

The  results,  expressed  in  terms  of  an  aerodynamic  radius  and  its  square 
which  is  proportional  to  lu*,*!  showed  that: 

a.  The  vertical  velocity  of  sedimentation  during  a time  has  an  inherent 
scatter  of  values; 

b.  The  average  aerodynamic  radius  of  a cylinder  is  somewhat  higher  than 
that  of  a unit  density  sphere  having  a radius  equal  to  the  equational 
axis  of  the  former; 

c.  The  standard  deviation  of  the  distributions  of  ^ and  ^ bounded 

by  the  values  appropriate  to  a flatewise  and  edgewise  fall,  is  between 
0.5^j  and  1 2p  of  the  average; 


d.  The  standard  deviation  for  cylinders,  for  example,  increases  with  the 
aspect  ratio  ( = length  over  diameter)  of  the  particle; 


i m 
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e.  The  average  velocity  of  sedimentation  of  cylinders,  for  example, 
does  not  change  with  , but  the  standard  deviation  decreases 
with  it. 

f.  The  "initial  history"  of  the  particle  is  forgotten  by  it  after 
times  of  sedimentation  of  the  order  of  a few  seconds. 

1 


II.  STATEMENT  OF  AIK  OF  REPORTED  STUDY 

In  view  of  the  existing  background  knowledge,  we  have  decided  to 
conduct  a research  along  the  following  lines: 

1 . Mobility 

a.  Continue  with  the  experimental  determination  of  the  principal 
resistance  coefficients  of  cylindrical  aerosol  particles  at  reduced 
pressures,  thus  increasing  the  upper  limit  of  the  Knudsen  number  of  our 
aerosol  system  and  treating  smaller  atmospheric  particles. 

b.  Investigate  experimentally  the  resistance  coefficient  of  cubical 
aerosol  particles  in  the  micron  size-range. 

2.  Deposition  in  still  air 

a.  Study  experimentally  deposition  in  still  air  of  cylindrical  aerosol 
particles,  and  compare  it  with  the  theoretical  xMonte  Carlo  calculations  of 
Gallily  and  Cohen  ( 1 976 ) . 


b.  Continue  the  Monte  Carlo  theoretical  calculations  (Gallily  and 
Cohen,  1976)  for  regularly-shaped  nonspherical  particles  so  as  to  find 
the  spatial  scatter  of  sedimentation  points  on  an  horizontal  plane. 


III.  EXPERIMENTAL 
Mobility 

A.  Cylindrical  Aerosol  Particles 

The  mobility  coefficients  of  cylindrical  aerosol  particles  at  a reduced 
pressure  were  experimentally  investigated  in  the  same  way  as  that  used  for 
their  determination  at  an  atmospheric  one  (Gallily,  1975). 

1 • Apparatus,  particles,  and  procedure 

The  apparatus  employed  for  this  purpose  was  essentially  the  same  as  that 
used  earlier  (Gallily,  1975)  but  for  the  incorporation  of  a vacuum  system 
attached  to  the  sedimentation  cell  and  glass  capillary  air  bleeders  for  the 
slow  introduction  of  the  aerosol  particles  to  be  determined  (Figs.  1 through 
4,  and  photograph  l). 


I I 

Fig.  1 . The  experimental  set-up  for  mobility  determination. 

F - stero-sedimentation  cell,  P - glass  fibers*  container, 

E - three  aerosol  introduction  bleeders,  A - vacuum  system  (pump, 
traps),  M - manometer,  D - stopcock 

[ 


} 
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Fig.  2.  The  stereo-sedimentation  cell  (vertical  cross-section) 

F - cylindrical  sedimentation  cavity,  G - scanning  electron  microscope 
stub,  , Bg-  viewing  microscopes'  port  holes,  E - Blocking, 


electromagnetically  activated  rod 


Fig.  3.  The  stero-sedimentation  cell  (horizontal  cross-section  at 
port-holes'  level). 


The  particles  were  glass  fibers  having  lengths  of  5-50  and  diameters  of 
0.3-1. 5 . These  were  produced  from  a Whatman  GF/C  Filter  Paper  raw 

material  (Whatman,  England)  by  an  electric  grinder,  and  dispersed  by  a 

shaking  action  as  described  earlier  (Gallily,  1975).  Their  density  at 

-3 

room  temperature  according  to  the  literature  was  2.23  g.cm  , which  was 
experimentally  checked  by  us  and  found  to  be  correct  within  -5?e. 

The  experimental  procedure  adopted  consisted  of  the  same  preliminary  steps 
of  system  levelling  and  calibration  as  used  in  the  atmospheric  pressure 
determinations,  evacuation  of  the  sedimentation  cell  F and  the  glass  fiber 
container  P (Fig.  1 ) till  a pressure  of  about  300  mm  Hg,  and  finally  a 
slow  controlled  introduction  into  the  cell  of  the  outside  air  which 
picked  up  fiber  particles  during  its  passage  through  P . The  trajectories 
of  the  particles  falling  into  the  lower  section  of  F (Fig.  2)  were  photo- 
graphed as  before  from  two  directions  of  sight,  the  location  of  the  sedi- 
mented particles  was  found  by  the  same  technique  reported  earlier  (Gallily, 
1975) > and  the  size  determination  was  performed  in  a scanning  electron 
microscope  (Stereoscan,  Mark  IV,  Cambridge,  England)  by  our  photogram- 
metric  method  (Gallily,  1975).  The  resistance  coefficients  pertaining  to 
motion  perpendicular  ( = K„  ) and  parallel  ( *»=*»)  to  the  cylinder's 

axis  were  calculated  according  to  the  mathematical  analysis  developed  by 


us  (Gallily,  1975). 


Here,  in  the  size  analysis  of  the  particles,  however,  they  were  not 
gold-shadowed  on  the  stubs  and  so  no  uncertainty  was  introduced  into 
the  measurement  of  their  binear  dimensions. 

2.  Results 

The  results  ere  given  in  Table  1 and  Fig.  5 in  the  form  of  the 
translational  resistance  coefficients  K,  (5  K,,)  and  £ = **0 
for  various  values  of  cylinder's  length  £ and  diameter  c at  reduced 
pressures.  is  the  resistance  coefficient  related  to  transverse  motion 
while  K33  is  the  coefficient  related  to  a longitudinal  one,  as  given  by  the 
diagonalized  translation  tensor 

K = ii  *«•  + K*1  + jfjj  C2,3 

/v 

the  theoretically  calculated  values,  kjj^are  based  on  the  continuous  fluid 
euations  (Happel  and  Brenner,  1965b) 

Arr 

K”  = c 
and 


2rr  ± c 


C 25  J 


where  = t/c  . 


o.s 

1 

CfM") 

H 

1/c 

K„(«XI0* 

(cm) 

(cm) 

***** 

(cm) 

(cm) 

(mm  Hg) 

t(°c) 

1 

58.2 

0.87 

67 

23.1 

17.5 

13.6 

8.3 

370.8 

20.3 

2 

24.7 

0.39 

63 

2.2 

1.5 

5.8 

3.6 

370.6 

21.7 

3 

11.1 

0.82 

14 

5.5 

4.1 

3.7 

2.5 

379.0 

22.8 

4 

17.0 

1.06 

16 

7.8 

6.6 

5.4 

3.6 

403.0 

21.2 

6 

44.6 

0.75 

60 

16.4 

9.6 

10.6 

6.5 

354.0 

21.6 

7 

8.9 

0.53 

17 

2.6 

1.6 

2.8 

1.8 

342.0 

25.0 

8 

25.1 

0.50 

50 

4.6 

3.8 

6.2 

3.8 

426.0 

26.2 

9 

19.7 

1.67 

12 

29.6 

20.8 

6.8 

4.6 

336.0 

28.3 

10 

5.9 

0.80 

7 

3.2 

2.1 

2.3 

1.7 

344.0 

28.3 

11 

46.2 

1.15 

40 

27.7 

19.4 

11.9 

7.5 

319.0 

24.9 

12 

21.3 

0.83 

26 

7.9 

6.6 

6.0 

3.9 

346.0 

28.1 

13 

40.4 

0.40 

101 

2.1 

1.2 

8.7 

5.3 

328.0 

26.7 

14 

1 14.8 

1.53 

10 

23.5 

14.2 

5.4 

3.8 

317.0 

29.3 

15 

} 13.7 

1.14 

12 

5.8 

4.9 

4.7 

3.2 

321.0 

28.9 

' 

16 

I 

j 22.9 

1.19 

19 

24.1 

22.0 

6.9 

4.6 

325.0 

29.5 

17 

118.9 

1.13 

17 

24.1 

11.8 

5.9 

3.9 

124.0 

29.3 

18 

5 58.3 

1 0.87 

67 

14.5 

9.3 

13.6 

8.3 

126.0 

29.3 

19 

1 38.2 

1 

j 0.73 

52 

10.9 

i 7.0 

! 

9.3 

5.7 

122.0 

28.0 

Table  1 : The  translational  resistance  coefficients  of  a cylinder  at  reduced 
pressures . 

K- . - experimentally  based  value;  X”-,; -t  - continuous  regime  theoretical 

value. 
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3.  Analysis  and  discussion 

The  experimentally  determined  resistance  coefficients  were  compared  to 
three  sets  of  calculated  values  (Table  2)  pertaining  to: 

a.  Continuous  fluid  coefficients, 

which  were  apriori  expected  to  be  higher  than  the  experimental  ones 
due  to  the  neglect  of  the  Knudsen  number  effect 

b.  Continuous  fluid  coefficients  for  short  cylinders  in  which 
correction  was  made  for  the  blunt  edges,  i.e.  finite  length,  of 
the  latter 

c.  Slip-continuum  regime  coefficients  for  a transverse  motion  of  an 
infinite  cylinder. 

The  corrected  resistance  coefficients  were  calculated  on  the  basis  of 
Heiss  and  Coull’s  study  (1952)  in  which  the  coefficients  were  expressed  in 
terms  of  the  empirically  determined  constant  G defined  by  Eq.  C 1 33  and  related 
to  K;i  by  Eq.  [l6]. 

Thus,  to  estimate  the  effect  of  the  blung  edge  of  a finite  cylinder  on 
its  translational  resistance,  one  can  use  Eq.  Cl 63  and  compare  the  coefficients 
of  the  cylinder  and  a prolate  ellipsoid  of  the  identical  dimensions 
according  to 


Kii  ej./k nt%K. 
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in  which  the  right  hand  bracketed  ratios  are  for  bodies  of  the 
axial  and  transverse  dimensions.  From  Heiaa  and  Coull's  study,  it  appears 
that  j Gt((  = 1 «06  for  l/c  = 5 and  that  it  approaches  1 for  a 

higher  length  to  diameter  ratio.  Thus, 


S «3./olv*0=  C 


2 r (c./zf  f 

(w/») 


'h 


which  is  the  actual  correction  factor  used  by  us. 

From  Tables  1 and  2,  it  is  seen  that: 

a.  The  translational  resistance  coefficients  of  the  cylinder  at  reduced 
pressures  ( K’-w  — 0.16)  are  of  the  same  order  of  magnitude  a8  those 
calculated  on  a continuous  medium  basis. 

b.  The  coefficients  for  cylinders  of  a diameter  above  1 ^*>1, 

approximately,  are  greater  than  the  calculated  ones  K',,*  , whereas  those 

/ 

related  to  cylinders  of  diameters  below  0,5  y'**  * approximately,  are  smaller 
than  , 

The  obvious  deduction  to  arise  was  that  the  effect  of  fluid  discontinuity 
(*^>0)  is  noticeable.  So  to  check  this  effect,  the  resistance  to  the 
transverse  translation  of  a cylinder  was  compared  with  the  theoretical 
equation  of  Pich  (1967)  in  which  the  Knudsen  number  influence  for  the 
slip-continuum  regime  has  been  taken  into  account 


Att  h-  

2. ooaz  - I.147  ^.^(1  + 0,7^^) 


C29J 


C27] 


- 24 


Table  2.  Comparison  of  resistance  coefficients  of  a cylinder  at  reduced  pressure*, 
Kiii*  “ experimental  values;  - continuous  fluid  values; 

\<u  c-  corrected,  calculated  values. 
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where  Fr,  is  the  resistance  per  unit  length,  and  ^ta_=  c/^. 

According  to  Eq.  [28], 

K„ 


A 7T  ^ 


2.00X0.  - fitj.-t-  l.ntin  K+i^-tr,  0 + O.n** 
for  a very  long,  essentially  infinite  cylinder. 

The  comparison  of  Eq.  [293  with  the  results  of  both  the  atmospheric  (Gallily, 
1975)  and  the  reduced  pressure  experiments  is  presented  in  Table  3*  from  which 
it  can  be  observed  that  for  cylinder's  diameters  above  0.5  j*'"'!,  Pich's  values 
are  lower  than  the  experimental  and  the  continuous  fluid  calculated  ones 
whereas  for  diameters  below  0. 5 »•»  the  former  approaches  the  experimental 
coefficients. 

To  smooth  out  the  fluctuations  in  the  determined  coefficients  , a 

regression  of  their  values  was  made  along  the  pattern  of  the  continuous  fluid 
expressions  (Eqs.[24]  and  [25])  as  the  theoretical  treatment  leading  to  those 
equations  was  expected  to  be  essentially  valid.  In  the  regression,  K';,-  values 
marked  in  Table  2 by  an  asterik  were  not  taken  into  account  according  to  a 
statistical  test. 

The  results  are  formulated  by 

Air 


K,.  = |.3I 


2*  + '/a 


3 


and 


Kn-  t.if? 


2tt 

Cn  s.4  - '/a 


C3IJ 


c(r-) 

X\\tP 

t* 

Hrt 

C C/o«) 

Kn/** 

Qjj 

IS 

5.45 

0.21 

0.59 

1.54 

0.82 

8.7 

0.49 

1.06 

2.69 

1.08 

6.1 

0.6 

0.76 

2.18 

3.64 

15.1 

0.74 

1.93 

4.51 

7.65 

15.6 

0.81 

2.02 

4.42 

9.98 

15.6 

0.63 

1.95 

4.45 

4.26 

86.4 

0.87 

11.28 

18.74 

74.22 

24.5 

0.83 

3.32 

6.95 

18.51 

56.2 

0.97 

7.44 

13*45 

**15.48 

0.72 

2.30 

2*87 

49.5 

1.42 

6.87 

13.12 

28.15 

0.99 

2.58 

1.12 

12.1 

0.76 

1.55 

3.83 

3.80 

16.5 

0.74 

2.12 

4.83 

5.22 

21.7 

0.96 

2.87 

6.34 

8.38 

58.2 

0.87 

7.48 

13.57 

23.09 

22.1 

0.70 

2.81 

5.98 

15.38 

24.7 

0.39 

2.82 

5.81 

2.17 

3.97 

1.0 

0.53 

1.94 

2.32 

11.1 

0.82 

1.42 

3.68 

5.46 

15.43 

0.71 

1.96 

4.54 

7.79 

17.1 

1.06 

2.26 

5.42 

7.85  i 

4.6 

0.75 

0.59 

1.93 

1.13 

8.3 

0.83 

1.06 

2.99 

2.92  j 

10.8 

0.82 

1.40 

3.60 

13.85 

44.6 

0.75 

5.60 

10.62 

16.44  1 

1 

13-8 

1.23 

1.88 

4.81 

10.44 

8.9 

0.53 

1.05 

2.78 

2.64 

| 20.3 

1.02 

2.70 

6.10 

13.57 

25.1 

0.5 

3.00 

6.18 

4.60 

[ 8.  7 

0.78 

1.12 

3.03 

3.23 

19.7 

1.67 

2.76 

6.76 

29.61 

5.8 

0.49 

0.70 

1.98 

0.73 

5.9 

0.8 

0.74 

2.31 

3.18 

12.2 

0.83 

1.59 

3.95 

2.75 

46.2 

1.15 

6.13 

11.88 

27.74 

18.8 

0.79 

2.43 

5.42 

4.95 

; 

21.3 

0.83 

2.71 

6.02 

7.87 

12.7 

1.07 

1.70 

4.35 

9.09 

, 40.4 

t 

0.4 

4.58 

8,75 

2.06 

58.3 

0.87 

7.45 

13.57 

14.55 

j 13.7 

1.14 

1.81 

4.67 

5.76 

j 

38.2 

0.73 

4.39 

9.34 

15.15 

22.9 

1.19 

3.06 

6.94  24.08  j 

1 

18.9 

1.13 

2.38 

5.93 

24.13 

Table  3.  A comparison  of  Pich’s  and  continuous  fluid  calculations  with  the 
experimental  resistance  coefficients  K (lj  w at  atmospheric  and 
reduced  pressures. 

- Pich's  values;  - continuous  fluid  value; 

3 

Kv,  t)<  - experimental  value,  all  multiplied  by  10  and  in  cm. 


Fig.  7.  The  regressed,  experimentally  determined  resistance  coefficient 
of  a cylinder  in  a parallel  motion  and  reduced  pressure. 

R - regressed;  T - continuous  regime. 

For  comparison's  sake,  the  equivalent  regressed  expressions  at  an  atmospheric 
pressure,  based  on  measurements  reported  earlier  (Gallily,  1975),  are  given 


10 
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Fig.  9.  The  regressed,  experimentally  determined  resistance  coefficient 
of  a cylinder  in  a parallel  motion  and  an  atmospheric  pressure, 

R - regressed;  T - continuous  regime. 

From  Eqs.  [30]  through'  [32]  and  Figs,  6-9  it  is  seen  again  that,  indeed, 
increasing  the  Knudsen's  number  (i.e.  reducing  the  pressure,  in  our  case) 
causes  a reduction  in  the  translational  resistance  coefficients,  as  expected. 

The  regressed  equations  can  obviously  serve  as  working  formulae  for  the 
introduction  of  the  experimentally  determined  tensor  K of  our  glass 
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fibers  into  the  appropriate  equation  of  motion  of  the  particles. 

As  a final  remark,  it  should  be  mentioned  that  for  the  cylinders  of 
diameters  above  1 whose  coefficients  were  found  to  be  higher 
than  the  continuous  fluid  value,  it  was  observed  that  none  of  them 
was  of  a perfect  shape;  furthermore,  it  was  noticed  that  the  less  perfect 
was  the  shape,  the  higher  was  the  value  of  K*;;  , which  leads  to  the  idea 
that  this  may  be  a factor  increasing  the  real  translational  resistance. 


B.  Cubical  Aerosol  Particles 

The  mobility  of  cubical  aerosol  particles  was  determined  in 
essentially  the  same  method  as  that  used  for  the  cylindrical  particles. 
Here,  however,  the  aerodynamic  behaviour  of  isolated  particles  in  a 
still  medium  is  simpler  than  the  behaviour  of  cylinders  since,  due  to 
the  symmetry  of  the  cubical  shape,  the  translation  tensor  K can  be 


written  as 


K I 

ry 


where  X is  the  idemfactor. 
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Following  Pettyjohn  et  al.  (1948),  one  can  express  u.*,^  by 
U*,*  = kj.  U-oo,$ 

where  U*>,s  is  the  Stokesian  velocity  of  a volumetrically  equivalent 
sphere  of  diameter  d*  . So,  denoting  the  side  of  the  cubical  particle  by 
Z,  one  gets  the  relationship 

ka.  = <8  23/  K cis 

in  which  is  hidden  through  K the  effect  of  the  fluid  discontinuity  on  the 
translational  resistance, 

1 • Apparatus,  particles. and  procedure 

The  apparatus  and  method  employed  for  the  mobility  determinations 
were  the  same  as  those  used  in  the  cylinders'  case. 

Though  falling  with  no  drift,  the  particles’  trajectories  were  photographed, 
as  beforehand,  from  two  directions  of  sight,  thus  facilitating  the  location 
of  each  measured  particle  on  the  stub  (Fig.  2)  and  increasing  the  accuracy 
of  the  velocity  determination. 

The  particles  themselves  consisted  of  sodium  chloride  (Frutarom,  Laboratory 
Chemicals,  Haifa)  cubes  produced  by  a generator  (Fig.  10)  which  was 
essentially  a spray-drying  instrument. 

In  this  instrument,  a 10^>  (w/w)  aqueous  solution  of  KaCl  was  dispersed  by 
a two-fluid  sprayer  A into  a 60  liter  plastic  drying  vessel  B.  The  sprayer 
was  activated  by  filtered  air  which  was  metered  by  F and  saturated  by  water 


Fig.  10.  The  spray-drying  generator. 

vapour  in  E so  that  the  nozzles  would  not  get  plugged  by  a crust  of 
dried  solid  (Gallily,  1969).  To  evaporate  the  produced  solution  droplets, 
another  stream  of  air,  metered  by  F and  dried  in  a SiO^  bed  D,  was 
introduced  into  B.  The  resulting  solid  particles  were  collected  in  the 
container  C and  used  for  the  mobility  determination.  These  particles  were 
nice  cubes  (Fig.  11 ) of  an  edge-dimension  of  0.5-2. 2 , which  was 

adjusted  by  choosing  an  appropriate  solution-concentration  and  spraying 
conditions. 
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The  density  of  the  cubes  at  room  temperature  was  taken  from  the  literature 
as  that  of  the  bulk  value  of  2.16(5)  g.cm  ,•  no  solution  filled  holes 
within  the  particles  were  assumed  because  of  crystallographic  energetical 
considerations  (Stdber,  1S75). 

The  linear  dimensions  of  three  edges  meeting  in  a comer  were  measured  for 
each  tested  particle  in  the  scanning  electron  microscope  mentioned  above, 
and  by  our  photogrametric  method  (Callily,  1975).  Only  particles  in  which 
any  edge  did  not  differ  in  length  from  the  other  two  by  more  than  lOfo,  were 
taken  into  consideration. 

The  dimensions  of  the  cubes  which  stood  this  "perfectness  test"  are  brought 
out  in  Table  4. 


V 

r 


Photograph  2.  A typical  sodium  chloride  cube  photographed  in  the  scanning 
electron  microscope.  Total  magnification  (on  the  photograph)  - 6460. 
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Concerning  electrical  charges  on  the  particles,  one  could  assume  that, 
due  to  the  spraying  mechanism,  they  did  exist  with  a distribution  which 
was  symmetrical  around  the  zero  charge;  however,  it  was  thought  that  these 
charges  did  not  influence  the  settling  velocity  of  the  particles  nor  did 
they  cause  an  increased  coagulation, as  could  be  perceived# 

Another  point  to  remark  was  that,  as  NaCl  particles  change  to  solution 
droplets  at  relative  humidities  of  about  70Je  in  room  temperatures,  they 
were  always  kept  in  a chemical  dessicator  over  a SiO^  bed. 


2.  Results 

The  translational  resistance  coefficients  * and  of  cubical 


aerosol  particles  are  given  in  Table  5 and  Figs.  11-12 


Table  4.  The  linear  dimensions  ( in  microns)  of  the  cubes  whose  mobility 


was  measured.  Z ” average  value  . 
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Table  5.  The  experimentally  measured  resistance  coefficients  of  translation 
of  cubical  aerosol  particles. 
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Fig.  11.  The  settling  velocity  of  cubical  particles  vs.  their  linear 
dimension. 
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3.  Analysis  and  dlscuasion 

To  smooth  out  the  fluctuations  of  the  experimental  velocity  determi- 
nations, a regression  was  made  on  the  basis  of 

fu*/  = 4 ^Z1  , C3aJ 

which  is  shown  as  a full  line  in  Fig.  1 1 . 

From  the  regressed  values  of  the  velocity,  one  could  obtain  a regressed 
expression  for  K (full  line  in  Fig.  12). 

The  coefficient  kA  calculated  according  to  Kq.  [37],  is  shown  in  Fig.  13 
together  with  the  regression  curve  (full  line)  which  is  based  on  the 
velocity  regression. 

The  curve  is  illustrative  in  comparing  the  relation  between  the  settling 
velocity  of  our  cubes  to  that  of  a volumetrically  equivalent  Stokesian 
sphere,  as  a function  of  2 . 

It  is  seen  that: 

a.  increases  with  2 

b.  For  Z 0.3  fX'mj is  approximately  2.8  times  greater  than 
ti.e  expected  values  of  Heiss  and  Coull  ( 1 952 ) at  continuous  fluid 
conditions,  but  only  1.4  times  greater  if  the  Cunningham  correction 
for  spheres  is  taken  into  account. 

c.  For  ZV  is  about  2.2  smaller  than  the  value  expected  from 
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the  above  study  of  Heiss  and  Coull 
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Fig.  13.  The  coefficient  to 
dimension  2 . 
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Here  one  should  note  that  this  behaviour  is  similar  to  the  behaviour 
found  in  our  study  for  cylinders. 

It  should  be  remarked  that  the  increased  resistance  to  motion  of  our 
cubes  over  the  valuee  found  in  liquids  in  earlier  studiee  can  not  be 
quantitatively  explained  by  assuming  a particle  density  lower  than  the 
bulk  one.  ho  fraction  of  solution  filled  holes  can  account  for  the  pheno- 
menon;  also,  the  external  shape  of  the  cubes  was  too  perfect  to  assume 
that  a density  decrease. 

One  should  note,  however,  that  our  cubes  were  not  perfect  bodies;  so,  the 
surface  defects  they  show  may  acccunt  for  an  increased  translational  resis- 
tance at  the  higher  end  of  the  Z-axia  (Fig.  13)  whereas  the  medium's  discon- 
tinuity (combined  with  the  surface  defects)  can  account  for  the  decrease  of 
K at  the  lower  end  of  Z . 

The  found  K valuee  can  be  introduced  into  the  equation  of  motion  L 1 j 
to  facilitate  the  solution  of  practical  problems  as  in  the  cylinder  case. 

Deposition 

The  deposition  of  nonspherical  particles  in  etill  air  can  be  measured 
sometimes  in  situ.  However,  in  many  inaccesible  systems,  such  as  the 
lower  parts  of  the  human  lungs,  one  has  to  infer  the  deposition  rate  from 
theory  (Gallily,  1975;  Gallily  and  Cohen,  1976)  once  the  eise  distribution 
of  the  particles  ie  known. 


The  sampling  of  particles  with  non-negligible  rate  of  sedimentation 
for  a size  distribution  analysis,  is  frequently  performed  by  utilizing 
their  settling  in  a quite  medium  on  an  horizontal  surface. 

In  the  case  of  spherical  particles,  one  can  relate  the  number  of  particles 
of  radius  r-  f f sedimenting  from  an  aerosol  atmosphere  in  time 
through  a clean  layer  of  thicmiess  AS  , to  their  number  per  unit  volume  of 
the  aerosol,  , by 

nr^  - ML  A+  - as) 


where  is  tile  appropriate  terminal  velocity  of  the  particles, 

M is  the  number  of  fields  of  view  checked,  and  L is  the  area  of  each 
field. 

ror  nonspherical  particles  the  relationship  is  more  complicated.  Thus, 
in  the  cylindrical  case,  Eq.  [ should  be  replaced  by 


~ WLNc.-c.CZ  C,  u U(  Ci  M -AS) 

' ' w.C,'  * / 

where  P © . u A+-  is  tile  probability  that  a cylinder  of  dimensions 

, */ 

L anc  c.;  will  have  in  time  At-  an  average  settling  velocity  u. 

c> 

This  prooability  has  been  calculated  by  Gallily  and  Cohen  (loc.cit). 

The  size  distribution  function,  X - N.  / hJ  is  given  now  by 

T({,e{  - / L.  c.  > 
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U%  r U;  -rr  o * 


/r  L~m  «W..  

(Tc.  ^ *.*/  Z Pt  (u;^-As)J 

i4.  = W*+  S c*  , • • 


In  our  experiments,  we  tried  the  technical  and  calculational  aspects  of 
a method  based  on  Eq. [41  3. 

I . Apparatus,  particles,  and  procedure 

The  apparatus  used  to  sample  the  particles  from  still  air  (Fig.  14) 
consisted  of  a sedimentation  ceU  having  a central  cylindrical  settling 
cavity  A,  a scanning  electron  microscope  collecting  stub,  and  two  lever- 
activated  diaphragm  stops  B of  0.1  mm  thickness,  one  at  a distance  of 
0.2  mm  from  the  glass-covered  stub  and  the  other  at  a distance  which 
could  be  varied  between  0.3  to  2.7  mm. 


Fig.  14.  The  sampling  cell  (scheoatical). 
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The  particles  were  glass  fiber  cylinders  of  the  same  origin  as  describee 
above  and  of  the  same  range  of  linear  dimensions.  These  particles  were 
produced  by  the  same  method  cited  beforehand  (Gallily,  1975).  The  glass 
covered  scanning  microscope  stubs,  however,  were  now  smeared  with  a thin 
layer  of  (filtered)  kerosene;  as  a result,  the  cylinders  which  sedimented 
were  pulled  down  by  surface  tension  action  towards  the  surface  and  lay 
flat  on  it,  as  previously  ascertained  also  for  similar  sampling  conditions 
by  Gallily  (1971 ). 

The  experiments  themselves  were  conducted  by  introducing  the  glass  fibers 
into  the  sampling  cell  with  the  aid  of  a Mariot  bottle  suction  while  the 
two  blocking  diaphragms  B (Fig.  14  ) were  shut,  opening  the  two  diaphragms 
at  time  t=0,  and  closing  the  lower  one  after  time  A.+ . 

The  sedimented  particles  were  measured  now  in  the  scanning  electron 
microscope  from  one  angle  of  sight  only,  which  saved  much  operation  time. 

2.  Regults 

The  times  of  sedimentation  were  always  26.6  seconds;  the  initial 
aerosol-free  layers  were  3.67,  2.66  and  0.33  mm,  giving  respectively 
minimal  average  velocity  of  sedimentation  C U.;  ^ iri  3 ^ 

of  138  100  (‘‘"’/sec  and  12.4  /*Vsec. 

The  dimensions  of  the  particles  found  on  the  stubs  under  the  conditions  of 
those  minimal  settling  velocities  are  given  in  Tables  6 through  8. 


Table  6.  The  linear  dimensions  of  cylindrical  particles  having 
a settling  velocity  equal  or  greater  than  138  pi*/  sec. 
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According  to  the  last  three  tables,  Eq.  l40J,  and  the  theoretical 
calculations  of  Gallily  and  Cohen  (1976)  for  the  probability 
appearing  in  Eqa.  [403-[4l],  one  gets  the  final  tables: 


Table  12.  ML  Ne.  for  various  values  of  and  CI/0  # 

u;  ^ ns  • 


10 

30 

50 

70  I 

0.5 
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_i 

Table  13.  M L N».  for  various  values  of  ?.>  and  Clo  . 

* I,  C,’  ' ' 
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Table  14.  W L N,  £ki  for  various  values  of  f;e  and  Cl(„ 
U(  ^ /2  . A /s«e  . 
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A a stated  above,  the  frequency  function  -f . can  be  now  calculated 

'i.c  i 

from  the  values  of  MlM.  ,,  . 

It  should  be  remarked  that  the  theoretical  calculations  of  P,  , 

w,  e,_  Uij&t 

were  baaed  on  continuous  fluid  values  of  the  resistance  coefficients. 

Thus,  for  the  case  of  Tables  12  ana  13,  no  particles  vithin  the  dimension- 
range  of  fi  = 10  - tOjr*  and  C;  = 0,5  - 0,5  could  settle  down  on  our 
stub  in  = 26.6  sec.  However,  one  should  remember  that  for  this  range 
of  sizes  the  resistance  coefficients  are  smaller  than  the  continuous 
medium  ones  and  so  this  should  have  to  be  taken  into  account. 

Judging  from  our  rough  results,  it  seems  that  the  bi-paraaetric  size 
distribution  function  -{ f . c of  the  glass  fiber  aerosol  has  a pronounced 
maximum  at  = 10  pin  and  C;(0  = 0.5  y* , and  that  it  slopes  down  towards 
larger  ti  and  C;  . 


IV.  THEORETICAL 

Spatial  distribution  of  sedimented  cylindrical  particles 

The  spatial  distribution  of  sedimented  cylindrical  particles 
on  an  horizontal  plane  was  calculated  by  a continuation  o f the  theoretical 
study  of  Gallily  and  Cohen  (Gallily,  1975;  Gallily  and  Cohen,  1976). 

Here,  however,  were  used  as  systems  of  coordinates: 
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An  internal  cartesian  system  whose  axes  are  parallel  to  the 
symmetry  axes  of  the  cylindr;  a minor  external  cartesian  system  in  which 
one  axis  is  parallel  to  the  direction  of  gravity,  the  other  lies  in  a plane 
determined  by  the  direction  of  gravity  and  the  polar  axis  of  the  particle, 
and  the  third  is  perpendicular  to  the  first  two;  a major  external  system  based 
on  the  horizontal  collecting  surface  and  in  which  one  axis  is  parallel  to  tne 
direction  of  gravity. 

1 . Results 

Typical  results  are  given  in  Figs.  15  through  16. 

The  calculations  were  performed  for: 

a.  Continuous  fluid  resistance  coefficients; 

b.  Continuous  fluid  resistance  coefficients  corrected  to  account  for  the 

1/3 

blunt  edge  by  the  multiplier  1.5  ; 

c.  Our  regressed  coefficients  in  an  atmospheric  pressure; 
a.  Slip-flow  coefficients,  given  by 

K„  = l<\,,e/0  + //O  C423 

and 

Ki»  = /(  I + a/c  ) t^3j 

where  X , tne  mean  free  path  of  the  air  molecules,  was  taken  to  be  O.Oo53  ffm. 

Observing  Figs.  Cl  53  through  Cl 83,  one  should  note  that  the  sedimenting 

particle  do  have  an  overall  drift,  as  found  in  an  earlier  study  (Gallily,  1974), 

and  that  this  drift  decreases  with  increasing  fall  distace. 


RADIAL  DISTANCE  XlO5 (cm.) 


60 


Mg.  16.  Cumulative  distribution  of  sedimented  cylindrical 
particles  vs.  the  radial  distance  from  the  origin  of  th« 
main  external  coordinate  system. 

Came  parameters  of  particles  and  calculations  as  given 
in  legena  to  rig.  15. 


Fie.  17.  Sedimentation  points  on  an  horizontal  plane  of  cylinders 

of  13.1  length  and  0.77  |*>v*  diameter  which  undergo  a rota- 
tional Brownian  motion.  Fall  distance  - 16.p5  mm.  Calculation 
parameters:  same  as  given  in  legend  to  Fig.  15. 

Average  aerodynamic  radius  - 1.37j/‘**i  ; standard  deviation  of 
aerodynamic  radii  - 0.07  p . 


' 25  3( 

RADIAL  DISTANCE  X I02  (cm.) 


Fig.  1b.  Cumulative  distribution  of  sedimented  cylindrical  particles 
vs.  radial  distance  from  the  origin  of  the  main  external 
coordinate  system. 


Same  parameters  of  particles  as  giver,  in  legei.d  to  Fig.  IT 
same  calculatioral  parameters  as  given  in  legend  to  Fig.  1 c 
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V.  GLOSSARY 


Cl,  - ellipsoidal  particle's  semi-axes 

a - radius  of  a spherical  particle 
A, -As-  - constants,  defined  in  text 

B - aerodynamic  mobility  of  a particle; 

B - mobility  tensor 
C - diameter  of  a cylindrical  particle 
C,  - number  concentration  of  aerosol  particles 
C^,  - factor,  defined  in  text 

ds  - diameter  of  a volumetrically  equivalent  particle; 

- cross  section  diameter,  defined  in  text 

D„  - rotational  diffusion  coefficient  of  an  axially  symmetric 

particle  about  an  equatorial  axis 

D,  O - operators,  defined  in  text 

JD/vt  - the  substantial  derivative  operator 

•{  - fraction  of  gas  molecules  diffusely  reflected  from  particle's 

surface  with  a temperature  accommodation  coefficient  of  1. 

4.  , - frequency  function  in  size  distriuution  of  cylindrical  particles. 

c;,c. 

F - force  ucting  on  a particle;  - external  force  ; F,  - 
fluiddynamic  force 
§ - gravitational  acceleration 

G - factor,  defined  in  text 

f,  j k - unit  vectors  along  particles'  principal  translation  axes 


1 
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Boltzmann's  constant 


k,  - coefficient  related  to  cubical  particles,  defined  in  text 
- Knudsen's  number,  defined  in  text;  k«n  - Knudsen's  number, 
defined  in  text 

K - translation  tensor;  - components  of  tensor  along  translation 
principal  axes  ( - coefficients  of  translational  resistance  )j 
K - coefficient  of  cubical  particles 
t - length  of  a cylindrical  particle 
L - area  of  an  observed  field  of  view 
Lc  - typical  length  of  a particle 


torque  acting  on  the  particle  with  respect  to  its  center  of 


mass  frl;  - external  torque;  L„  F - fluiddynamic  torque 

inn  - mass  of  an  aerosol  particle 
yn  - mass  of  fluid  displaced  by  a particle 

M - number  of  field  of  views  checked  in  size  distribution  analysis 
- number  of  sedimented  spherical  particles  of  radius  rj 
which  are  counted 

Nr  - number  density  of  spherical  aerosol  particles 


'Vef 


number  of  sedimented  cylindrical  aerosol  particles  which  are 
counted 


Nl,.  r.  - number  density  of  cylindrical  aerosol  particles 
p - total  pressure  in  gaseous  medium 
P - stress  triadic  defining  B 


WWW, 


■UWft'WU? 
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"3* 


(°-y,(X)  - Legendre's  polynomial  of  degree 


P - probability  of  a cylindrical  aerosol  to  have  with  the 

'#  * 


time  let  an  average  settling  velocity  u«.  , normalized 

W | c ^ 

U,  a ' 


SO 


that  £ P =1 

Ula4V«tt*'Ci  •" 

Qi  - particle's  instrinsic  diadic,  defined  in  text 

- radius  vector  from  center  of  mass  to  a surface  point 
R - radius  of  adjusted  sphere 

(?tt  - translational  Reynolds  number,  defined  in  text;  R&  r - 
rotational  number,  defined  in  text 
Rs  - particle' 8 Stokes  radius;  f?a-  par  tide's  aerodynamic  radius 
JS  - surface  element  of  a particle 
£ - surface  of  a particle 

AS  - thickness  of  aerosol  - free  layer 
■fc  4_  - dummy  variable 
~t  - time;A+  - time  difference 

T - temperature  (absolute  scale),  - surface  temperature 

- translational  velocity  of  particle's  center  of  mass;  - 

terminal  sedimentation  velocity  ( Uo»(y  > U*>,i  " components 

of  velocity  along  external  axes);  u.  (t  U;)  - sedimentation 

w.c,* 

velocity  of  a cylindrical  particle 
V - fluid  velocity 

/v 

- cartesian  coordinates  along  particle's  translational  axes; 

X, $,i  - external  coordinates  (i  - parallel  to  direction  of 

gravity) 


X - cos  6 ($  - colatitude  angle  between  the  polar  axis  of  an 

axially  symmetric  particle  and  an  initial  direction  in  space 
Z - length  of  cube's  side  ( three  sides  meeting  in  a 

comer)  ;Z  - average  value. 


Greek  Letters 


r*  - 

<p\ar’  - 

A - 
A - 
> - 

Y-  - 
1 - 

- 

IT  _ 
« 

Sf  - 
- 

4>  = 
t - 
co  - 

/■v 


integrals,  defined  in  text 

direction  cosines  with  respect  to  particle's  translation 
defined  in  text 
dummy  variable 

mean  free  path  of  gaseous  molecules 

dynamic  viscosity  of  fluid 

kinematic  viscosity  of  the  fluid 

integral,  defined  in  text 

stress  tensor  within  a fluid 

density  of  a particle 

tangential  stress  on  particle's  surface 

l/c 

sphericity  factor 

rotational  velocity  of  the  particle 


WHUiiljyull  llliliiili_l  WF>  U.IIIHPI-III..I 
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Subscripts 

C - continuous  regime  value 
e - external 

F - fluiddynamic 

i - i-th 

M - center  of  mass 

o - central  value 

t - tangential 
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